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Abstract. We have performed a first fully 3-D GRMHD simulation 
with Schwarzschild black hole with a free falling corona. The initial sim- 
ulation results show that a jet is created as in previous axisymmetric 
simulations. However, the time to generate the jet is longer than in the 
2-D simulations. We expect that due to the additional azimuthal dimen- 
sion the dynamics of jet formation can be modified. 



Introduction 



Relativistic jets have been observed in active galactic nuclei (AGNs) and micro- 
quasars in our Galaxy, and it is believed that they originate from the regions 
near accreting black holes. To investigate the dynamics of accretion disks and 
the associated jet formation, we use our newly developed full 3-D GRMHD code. 

Recently, Koide, Shibata, & Kudoh (1999) have investigated the dynamics 
of an accretion disk initially threaded by a uniform poloidal magnetic field in 
a non-rotating corona (either in a state of steady fall or in hydrostatic equilib- 
rium) around a non-rotating black hole. The numerical results show that the 
disk loses angular momentum by magnetic braking, then falls into the black hole. 
The infalling motion of the disk, which is faster than in the non-relativistic case 
because of the general-relativistic effect below 3rg (rg is the Schwarzschild ra- 
dius), is strongly decelerated at the shock formed by the centrifugal force around 
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r = 2rg by the rotation of the disk. Plasmas near the shock are accelerated by 
the J x B force, which forms bipolar relativistic jets. Inside this magnetically 
driven jet, the gradient of gas pressure also generates a jet over the shock re- 
gion {gas-pressure-driven jet). This two-layered jet structure is formed both in 
a hydrostatic corona and in a steady-state falling corona. Koide et al. (2000) 
have also developed a new GRMHD code in Kerr geometry and have found that, 
with a rapidly rotating (a = 0.95) black-hole magnetosphere, the maximum ve- 
locity of the jet is 0.9 c and its terminal velocity 0.85 c. All of the previous 2-D 
GRMHD simulations described here were made assuming axisymmetry with re- 
spect to the z-axis and mirror symmetry with respect to the plane z = 0; the 
axisymmetric assumption suppressed the azimuthal instabilities. 

2. 3-D GRMHD Simulations: Equations and Numerical Techniques 

Our basic equations are those of Maxwell for the fields and a set of general- 
relativistic equations representing the plasma, namely the equations of conser- 
vation of mass, momentum, and energy for a single-component conducting fluid 
(Weinberg 1972; Thorne et al. 1986). In making the simulations, we use these 
equations with the 3+1 formalism (for details, see Koide, Shibata, & Kudoh 



3. Preliminary 3-D GRMHD simulations with a Schwarzschild black 



In order to investigate how accretion disks near black holes evolve under the 
influence of accretion instabilities such as the magnetorotational instability, the 
use of a fully 3-D GRMHD is essential. 

3.1. Initial and boundary conditions 

In the assumed initial state, the simulation region is divided into two parts: a 
background corona around a black hole, and an accretion disk (Fig. la and lb). 
The coronal plasma is set in a state of transonic free-fall flow, as in the case of 
the transonic flows with T = 5/3 and H = 1.3; here the sonic point is located 
at r = 1.6rg. The Keplerian disk in the corona is set in the following way. The 
disk region is located at r > r-Q = 3rg, |cos6>| < 5 = 1/8. Here the density is 
100 times that of the background corona (Fig. la), while the orbital velocity 
is relativistic and purely azimuthal: = vk = c/[2(r/rs — l)] 1 ^ 2 - (Note that 
this equation reduces to the Newtonian Keplerian velocity v<p = \JGMjr in the 
non-relativistic limit rg/r <C 1). The pressure of both the corona and the disk 
are assumed equal to that of the transonic solution. The initial conditions for 
the entire plasma around the black hole are: 
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, x J (0,0, wk) (r>r D and |cot(9| < <5) ( ) 

(Vr,ve,Vi)-\ { _ VScA0) (r<r D or |cot0| > 6) {6) 

where we set S = 0.125; the smoothing length is 0.3rs- 

In addition, there is a magnetic field crossing the accretion disk perpendicu- 
larly. We set it to the Wald solution (Wald 1974), which represents the uniform 
magnetic field around a Kerr black hole: B r = Bqcos9, Bq = —aBosinO (where a 
is the lapse function, a = (1 — rs/r) 1 / 2 ). At the inner edge of the accretion disk, 
the proper Alfven velocity is va = 0.015c in a typical case with Bq = 0.3\/ poc 2 , 
where the Alfven velocity in the fiducial observer 

v A = B/^p+[Tp/(T-l) + B 2 ]/c 2 . (4) 

The plasma beta of the corona at r = 3rs is [3 = p/B 2 = 1.40. The simulation 
is performed in the region l.lrs < r < 20rs, < 9 < tt, < (/) < 2ir with 
100 x 120 x 60 meshes. The effective linear mesh widths at r = l.lrs an d at r = 
20rs are 5.38 x 10 _3 rs and 0.97rs, respectively, while the angular spacings along 
the polar and azimuthal directions are 5.2 x 10~ 2 rad. A radiative boundary 
condition is imposed at r = l.lrs and at r = 20rs: 

= u™ + < +1 - n™, (5) 

where the superscripts n + 1 and n denote the time steps and the subscripts 
and 1 refer to the boundary and to its neighbor meshes, respectively. The 
computations were made on an ORIGIN 2000 computer with 0.898 GB internal 
memory, and they used about 47 hours of CPU time for 10000 time steps with 
100 x 120 x 60 meshes. 



3.2. Simulation results 

Figure 1 shows the evolution of 3-D simulation performed in the region l.lrs < 
r < 20r s , < 9 < tt, and < <j> < 2vr with 100 x 120 x 60 meshes. The 
parameters used in this simulation are the same as those of the axisymmetric 
simulations shown in Fig. 6 of Koide, Shibata, & Kudoh (1999). In this figure, 
the colored shading shows the proper mass density ((a) and (c)) and the pressure 
((b) and (d)) on logarithmic scales; the vector plots show the magnetic field ((a) 
and (c)) and flow velocity ((b) and (d)). 

The black circle represents the black hole. Figures la and lb present the 
initial conditions, which are the same as in the 2-D simulations (Koide, Shibata, 
& Kudoh 1999). At t = 39.2r s comparing with Fig. 6c (t = 40.0r s ) in Koide 
et al. (1999) the jet is less generated. At t = 60.0rs the jet is generated as 
2-D simulation at the earlier time {t = 40.0rs). At t = 73.9rs, the jet is clearly 
created around r = 4.5rs, which is shown by the enhanced density (Fig. lc) 
and pressure (Fig. Id). As in the 2-D simulations the jet is generated in a 
hollowed cylindrical form. The magnetic field is twisted by the accretion disk 
and pinched, which increases the magnetic field pressure and generates the jet 
near the black hole. The delay of jet formation seems to be due to reduction of 
shock formation at r = 2rs caused by the additional freedom in the azimuthal 
direction. Further investigation will be reported elsewhere. 
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Figure 1. For the fully 3-D simulation, these panels present the time 
evolution of the proper mass density with the magnetic field (B X ,B Z ) 
((a) and (c)) and the proper pressure with the flow velocity (v X iVz) 
((b) and (d)) in a transonic free-fall (steady-state falling) corona with 
an initially uniform magnetic field, at t = O.Ors ((a) and (b)), and 
t = 73.9rg ((c) and (d)). The jet is formed around r = 4.5rs as in the 
2-D simulation. 

4. Discussion 

This simulation result is initial and we will perform more simulations and in- 
vestigate effects of the third dimension. Further results will be reported else- 
where. K.N. is partially supported by NSF ATM 9730230, ATM-9870072, ATM- 
0100997, and INT-9981508. The simulations have been performed on ORIGIN 
2000 at NCSA which is supported by NSF. 
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